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Abstract. A cosmological model with anisotropic dark energy is analyzed. The amount of 
deviation from isotropy of the equation of state of dark energy, the skewness 5, generates an 
anisotropization of the large-scale geometry of the Universe, quantifiable by means of the actual 
shear Eo. Requiring that the level of cosmic anisotropization at the time of decoupling is such to 
solve the "quadrupole problem" of cosmic microwave background radiation, we find that \8\ ~ 10~ 4 
and | So | ~ 10 , compatible with existing limits derived from the magnitude-redshift data on type 
la supernovae. 
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I. INTRODUCTION 

The isotropy of the cosmic microwave background (CMB) radiation, first seen by the Cosmic Background Explorer 
(COBE) satellite [l| and then reinforced by the Wilkinson Microwave Anisotropy Probe (WMAP) data 0, together 
with the assumption that we are not in any special position in the Universe, underlies the Cosmological Principle, 
according to which we live in a homogeneous and isotropic universe described by a Robertson- Walker line element. 

Tiny deviations from perfect isotropy, at the level of 10~ 5 , have been also reported by COBE [3| and thereafter 
confirmed by the high resolution WMAP data. The observed CMB anisotropy spectrum is in impressive agreement 
with the predictions of the ACDM concordance model, namely the widely accepted standard model of (inflationary) 
cosmology [4|. 

Nevertheless, data from 7-yr WMAP observations [ff do not have removed an effect seen at lower significance in 
the COBE data, that is a lack of power on the CMB anisotropy quadrupole moment. 

Since measurements at the largest angular scales are affected by foreground and other systematic effects, we have to 
wait for the full results of the PLANCK mission Q that will check, using a different frequency coverage and scanning 
strategy than WMAP, the anomalously low quadrupole so far observed. 

Surprisingly, however, analyzing the first light sky map released by the PLANCK team, which covers about only 
the 10% of the sky, Liu and Li [7[ have recently claimed that the amplitude of CMB power spectrum on large scales 
is significantly lower than that reported by the WMAP team. 

The "smallness" of CMB quadrupole amplitude, referred to as quadrupole problem, signals a suppression of power 
at large scales and has been intensively studied |8| meanly because it may signal either a nontrivial topology or a 
deviation from isotropy of the large scale geometry of the Universe. Indeed, recently enough, it has been shown that 
a spatially homogeneous but anisotropic cosmological model of Bianchi type I described by a plane-symmetric line 
element, the ellipsoidal universe 0, allows a better matching of the large-scale CMB anisotropy data. 

Various mechanisms could have induced a planar symmetry in the spatial geometry of the universe: topological 
defects (e.g. cosmic stings, domain walls) [HI], a uniform cosmological magnetic field [fjlicj, magnetic fields possessing 
planar symmetry at cosmic scales , or a moving [3] dark energy. 



In particular, the mechanism proposed by Koivisto and Mota of dark energy with anisotropic equation of state 14 



HI (see also Refs. ^$2$) is very attractive, since cosmic anisotropy originates from the actual dominant component 
of the Universe and then could be directly tested, for example, by either observations of magnitude and redshift of 
type la supernovae or cosmic parallax effects of distant sources [211. 

The aim of this paper is to connect the dark energy anisotropization mechanism to the ellipsoidal universe proposal. 
As we will see, the deviation from isotropy of the equation of state of dark energy triggers an anisotropization of the 
Universe which can be well described by an ellipsoidal universe model. The level of such an anisotropization can be 
then chosen to solve the quadrupole problem. 

The plan of the paper is as follows. In section II, we briefly review the ellipsoidal universe proposal. In section 
III, we introduce a cosmological model with an anisotropic dark energy component and we study its connection with 
both the ellipsoidal universe model and the quadrupole problem. Finally, we draw our conclusions in section IV. 
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II. ELLIPSOIDAL UNIVERSE: CMB QUADRUPOLE 

Planar symmetry, in a cosmological context, is described by the Taub line element [22j ]: 

ds 2 = di 2 - a 2 {t){dx 2 + dy 2 )- b 2 (t) dz 2 , (1) 

where the scale factors a and b are normalized as a(to) — b(to) — 1 at the present (cosmic) time to. 

As shown in Ref. |9(, the above metric generates a quadrupole term in the CMB radiation, through the Sachs- Wolfe 
effect, which adds to that caused by the inflation-produced gravitational potential. 

More precisely, let us introduce the CMB temperature anisotropy in terms of spherical harmonics Y/ m , 

oo / 

AT(n) = a 'mYi m (n), (2) 

1=2 m=—l 

where n is a unit direction vector. The coefficients of the expansion, a; m , are related to the angular power spectrum, 
Ci, through 

7U— — 1 

and the quadrupole (I = 2) component is defined as 



Q=^C 2 . (4) 
The observed CMB anisotropy map is then a linear superposition of two contributions [23j]: 1 

O-lm — a^n, + a \mi (5) 

where af^ comes from the anisotropic spacetime background, while the a\ m term is the standard isotropic fluctuation 
due to inflation-produced gravitational potential at the last scattering surface. 
The 7- year WMAP data give, for the quadrupole amplitude, the value Q 

Q2 bs ~ 200^K 2 , (6) 
while their best-fit for the isotropic ACDM standard cosmological model gives 

Qi~f200^K 2 . (7) 
Even taking into account the cosmic variance [2~ij ]. which dominates theoretical uncertainties at low I, 

|q^~ 759 M K^, (8) 

the quadrupole amplitude remains anomalously low. 

The quadrupole anisotropy associated to the anisotropic spacetime background has been instead calculated in 
Ref. 0: 

2 

2a = — ^ ^cmbedec , (9) 

where T cm b ~ 2.73 K is the actual (average) CMB temperature and edec = e(^dec) is the "eccentricity" 



V 7 !- (b/a) 2 , a>b, 

^l-(a/b)\ a<b, [U> 



As pointed out in Ref. Il5ll . this is an approximation which is good when the spacetime background anisotropy is small (this will be 
indeed our case, see section III), and that a rigorous treatment of CMB anisotropics should presuppose the analysis of perturbations in 
the anisotropic background Eq. |[T}. 
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evaluated at the time of decoupling idee- 

The "ellipsoidal universe" proposal could indeed explains the "low" value of the observed quadrupole, Eq. ([6]), 
compared to that predicted by the standard cosmological model, Eq. ([7]). This, however, requires a suitable orientation 
of the anisotropy associated to planar geometry with respect to the inflation-produced one, in such a way to lower 
the overall power to a sufficient extent. As shown in Refs. 13, ll| this is attained for eccentricities approximatively 
given by 



a5[3V73-5sgn(q-6)] Q 1 
Cdcc " 24 T cmb ' [ii > 

where sgnx is the sign function (sgnx = ±1 if x ^ 0). In the Appendix we show that the probability that such an 
orientation occurs by chance is, indeed, not negligible. 

In the next Section, we show that an anisotropic dark energy component causes an anisotropization of the Universe 
described by the line element (TTJ. In particular, the skewness, which parameterizes the deviation from isotropy of the 
equation of state of dark energy, and the shear, which quantify the level of cosmic anisotropization, will be connected 
to the amount of eccentricity at decoupling. 

III. ANISOTROPIC DARK ENERGY: SKEWNESS AND SHEAR 

The most general energy-momentum tensor consistent with planar symmetry is 

T"„ = diag(p,-pll,-pH ! -p x ), (12) 

where p is the energy density, while p" and p^ are the "longitudinal" and "normal" pressures. Taking into account 
the energy-momentum tensor (|12p . the Einstein's equations read 

3(1 - S 2 )i? 2 = 8ttGp, (13) 
3(1 - £ + T, 2 )H 2 + [(2 - E)£T]- = -8ttG P \ (14) 
3(1 + T,) 2 H 2 + 2[(1 + £)#]' = -8<jrGp x , (15) 

where a dot denotes a differentiation with respect to the cosmic time. Here, we have introduced, in the usual way, 
the cosmic shear, S, and the "mean Hubble parameter", H, as 

^ = {H a -H)/H, H = A/A, (16) 

where H a = a/a and A = (a 2 b) 1 ^ 3 is the "mean expansion parameter". 

Instead of considering Eqs. (fl"5)) and (|T4|) . we will analyze the equation obtained by subtracting side by side them, 

3(i?S)' +9# 2 £ = 8ttG(p 11 -p- 1 ) (17) 

(from which it is apparent that the source of the shear is proportional to the difference between the longitudinal and 
normal pressures of the anisotropic fluid), and the equation coming from the conservation of the energy- momentum 
tensor, T^ u .^ = 0, which gives 

P + 3g (p + 2P " 3 P± ) + 2ff (p'l -p ± )Z= 0. (18) 

To proceed further, we assume that the anisotropic fluid defined by the energy-momentum tensor (|12j) is indeed made 
up of three different components: an isotropic radiation component (r), an isotropic dark matter component (m), and 
an anisotropic dark energy (DE) component, 

P = pr + Pm + PDE, (19) 
p 11 =p r +p m +pJjj E , (20) 
V 1 " =Pr +Pm+PDE) (21) 

with equations of state: p r = p r /3, p rn = 0, and 

Pd E = w^de, Pde = w^Pde- (22) 
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Moreover, we will assume that the w" and coefficients are constants and that the three components are noninter- 
acting. The latter assumption ensures that each component is separately conserved, so that Eq. (|18p gives 

Pr =p^A-\ Pm = pV>A- 3 , (23) 
Pde + [3(1 + w) + H PDE = 0, (24) 

where from now on an index "0" defines quantities evaluated at the actual time, and we have introduced the "mean 
equation of state parameter" w and "skewness" S as 

w = (2W 11 +w L )/3, S^w^-ur 1 . (25) 

Finally, we assume that S and S are small quantities (as we will verify a posteriori) so that we can neglect the second 
term in the square brackets of Eq. (|24|) . In this case, Eq. (f24|) simply gives: 

p DE =p^- 3(1+,u) - (26) 

Hereafter, we use this approximate result for the evolution of dark energy density. 
Introducing the "mean density parameters" 

where X = r,m, DE, and taking into account Eqs. ([!"§)) - (|2"T1) . the shear equation (fP7)) gives 
where 



H(A)/H = y/n r A-* + Q m A~ 3 + n DE A-^+^ (29) 

and we have defined the function 

dx 



E(A) = DE / ~r 
Jo % 



H(x)/H Q 



(30) 



It is worth noting that the density parameters are not all independent, sice evaluating Eq. (|29[) at the present time 
gives 

O r + O m + DE = 1- (31) 

In order to not to spoil the predictions of the standard isotropic model (such as those coming from Big Bang Nucle- 
osynthesis and Large Scale Structure formation), we only consider anisotropic cosmological models which isotropize 
at early times. Therefore, we impose the isotropization condition: 

lim = 0. (32) 

Taking into account Eqs. (|2"51) - (|3"0")) . it is easy to verify that the above condition is satisfied if and only if 

So = £cA (33) 

where we have taken into account that the dark energy component is subdominant with respect to the radiation 
one for A — >• 0. For completeness, we give the asymptotic expansion of the shear for small values of the expansion 
parameter: 

A«l: EW^^A 1 "*-. (34) 
Moreover, since dark energy dominates over the other components in the far future, A — > oo, we get from Eq. (|28[) : 

= lim = 26 , (35) 

A-t-oo 3(1 — W) 
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where we used Eq. ([551) . 

Finally, taking into account Eqs. 
holds: 



(|30p and Eqs. ((55)) - (1551) . it is possible to show that the following inequality 



max |E0A)|=|E c 

Ae(0,oo) 



(36) 



Equations ()33[) -(j36 [) makes evident the fact that an anisotropy in the equation of state of dark energy (5) generates an 
anisotropy in the cosmic geometry (E), and that the Universe will never isotropize, although its level of anisotropization 
is very low being S a small quantity. 

In Fig. 1, we plot the shear as a function of the expansion parameter for various values of w and Ode- Here and in 
the following we take O r = 0.83 x 10" 4 [25| . 

To proceed further, we note that the eccentricity and the shear are connected by the following relation: 



f A dx 

6sgn(a — b) / — E(x), 
Ji x 



(37) 



valid for small eccentricities, e«l, and coming from definitions (|T0f and (|T6)) . The above equation calculated at the 
time of decoupling, together with Eqs. ([2"8]) and J55J), implies that sgn(a — b) = — sgn<5 and allow us to write S and E 
in the form: 



\6\ = c s (w,n BE ) e 2 dcc , 
|E | = c s (w,^de) e\ ec , 



where 



eg 



r 1 dx E(x) 1 

J AdDc x* H(x)/H _ 

E Q C S . 



(38) 
(39) 



(40) 
(41) 



Here, ^4dcc = ^.(^dcc) is the mean expansion parameter evaluated at the time of decoupling. In the following we simply 
assume that Ad cc = 1/(1 + Zd cc ), where Zd oc — 1090 [3] is the redshift at decoupling. 

In Fig. 2, we plot the skewness and the actual shear as a function of the mean dark energy density, for various 
values of the mean dark energy equation of state parameter. It is evident from the figure that, for a wide range of 
values of Ode and w, the parameters S and Eq are indeed small quantities, as we have previously assumed. 
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FIG. 2: The skewness S (upper panel) and the shear So (lower panel) as a function of the mean dark energy density S7de for 
different values of the mean dark energy equation of state parameter w. The values of 5 and Eq ensure that the eccentricity at 
decoupling is just that given by Eqs. (0 and 



The smallness of the parameters of anisotropy, S and So, is a reason to believe that the values of mean parameters 
firjE and w are very close to the analogous ones for the isotropic standard cosmological model, which according the 
7-yr WMAP data are [Hj]: 

^(isotropic) = Q ± 029 ^ g8% c L ^ ^ 

^(isotropic) = _ .980 ± 0.053 (68% C.L.). (43) 
From the above equations and looking at Fig. 2, we finally get an order of magnitude estimate of 5 and So : 

|<5|~10- 4 , |£o|~10- 5 . (44) 



The above value of 5 is compatible with the constraint obtained by Koivisto and Mota [14j coming from the analysis 
of luminosity distance-redshift relation of type la supernovae. 2 Instead, no constraint exists on Eq in the literature. 



2 The dark energy equation of state parameter w and the skewness <5 introduced in this paper correspond, respectively, to the parameters 
w + 7 and -37 defined in Ref. [Jj|. Translating the results of Ref. 14] to our case, we get -0.15 < <5 < 0.21 (68.3% C.L.). 
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Nevertheless, a full analysis of magnitude-redshift data on type la supernovae can put limits both on S and So and, 
indeed, preliminary results [27j indicate that 

-0.40 < 5 < 0.13 (la C.L.), (45) 
-0.026 < E < 0.014 (la C.L.). (46) 

This leaves open the possibility of having an ellipsoidal universe via an anisotropic dark energy. 



IV. CONCLUSIONS 

The 7-year WMAP observations H do not have alleviated the so-called quadrupole problem of CMB anisotropy 
spectrum, and a tension between data and the quadrupole amplitude predicted by the best-fit ACDM concordance 
model still remains. 

On the other hand, an anisotropic cosmological model described by a Bianchi type I metric, the "ellipsoidal uni- 
verse" [f| [l(| , allows a lower quadrupole and better matches the large-scale CMB anisotropy data. 



In the seminal papers [lj, |15[ by Koivisto and Mota, the peculiar features of a universe filled with an anisotropic 
dark energy fluid were deeply and exhaustively studied. 

In this paper, instead, we have investigated the possibility that a dark energy component with anisotropic equation 
of state could generate an ellipsoidal universe with the right characteristics to explain the low quadrupole in the CMB 
fluctuations. 

The amount of anisotropy in the equation of state of dark energy, the skewness 5, is transferred to the background 
geometry which gets anisotropized at a level that today is given the actual shear E . Imposing that such a level of 
cosmic anisotropy is exactly that necessary to generate a quadrupole term in the CMB spectrum in such a way to 
match the "low" value of the observed one, we have obtained the estimates: |<5| ~ 10~ 4 and |SqI ~ 10~ 5 . 



These values are well within the la confidence region allowed by magnitude-redshift data of type la supernovae [14 . 
|27| |. whose analysis constitutes, up to today, the only test available in the literature constraining cosmic anisotropy 
of Bianchi type I. 

This kind of test when combined with other tests on cosmic anisotropy, such as that coming from the study of CMB 
polarization spectrum, could put severe limits on the existence of anisotropic dark energy, confirm the ellipsoidal 
universe proposal or rules it out. All of this is, however, beyond the aim of this paper and will be the subject of future 
investigations. 



APPENDIX 

The total CMB quadrupole intensity, Qtot, is given by Eqs. ©-([5]) as [{j Hoj: 

QL = Ql + Q 2 1 -2fQ A Q 1 , (47) 

where Q\ and Qa are respectively given by Eqs. (O and (|9]), 

1 + 3 cos(2tf) + 2^6 sin$ [sin-i?cos(2(^ + ^3) — 2cos"i9cos((p + ^2)] , . Q \ 

j(y,tp<p2,<P3) = — y= , (48) 

and we are considering, for the sake of simplicity, the case a > b. In Eq. (J48J , $ £ [0, tt] and tp £ [0, 2ir] are the 
polar angles defining, in the galactic coordinate system, the direction of the axis of symmetry associated to planar 
symmetry of ellipsoidal universe model foL lioj. Together with the eccentricity at decoupling, edoc, they completely 
define the coefficients a^m- The parameters (^2,3 £ [0,27r], instead, are unknown phases (which, roughly speaking, 
define the "direction" of the inflation-produced quadrupole) that together with Q\ completely define the coefficients 
a\ m . Explicitly, we have P. Hoj: 

40 = -^[l + 3cos(2tf)]e;L c , 

41 = = sin(2tf) 4 C , (49) 

4 2 = (4,- 2 )* = -\/^e- 2 ^sin 2 tfei cc , 
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FIG. 3: The probability distribution function of / (upper panel), and the probability that the observed quadrupole Q obs is in 
the la confidence interval of the theoretical total quadrupole Qtot i a t the varying of the eccentricity at decoupling e doc (lower 
panel) . 



and 

a 2o = d 3 2i , 

4i = -(4,-i)* = y|e i * a Ci ) (50) 

Let us assume that there exists some mechanism able to produce an anisotropization of the Universe, so that edcc is a 
parameter fixed by the model. Generally, on the other hand, the direction of the symmetry axis {§, ip) cannot be fixed 
by the model and should be considered as unknown. The same its true for the direction of the inflation-produced 
quadrupole (<f>2,<f>3)- We therefore assume that z9, ip, (j>2, and </>3 are stochastic variables which assume random values 
in their intervals of existence. In Fig. 3 (see the upper panel), we show the probability distribution function of /, 
F[f], derived by performing a Monte Carlo simulation. 

We want now to calculate the probability, whatever are the values of i?, ip, 02, and </> 3 , with ed cc being fixed, that 
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the observed quadrupole Q 2 hs , given in Eq. (j6)), is in the la confidence interval of Qf ot , namely 



P(e dec ) = P(e dcc | Q 2 ohs e [Q 2 tot - <r tot , Q 2 tot + a tot ]). 



(51) 



The uncertainty <7tot on Q 2 is obtained by simply propagating the uncertainty on the inflation-produced quadrupole 
Qi, i.e. the cosmic variance in Eq. (jH)), by means of Eq. (|47|) : 



O"tot 



Qa 
J Qi 



(52) 



It is easy to see that the condition Q 2 hs € [Q 



2 

tot 



f± (edec 



ctot , Qtot + "tot] is equivalent to the condition / e [/_, / + ], where 
1 ± 



si 



^obs 



2T 



QaQi 



We then have that 



-P(edec) 



/+ 



dfF[f] 



(53) 



(54) 



f- 



gives the probability, at fixed edec, that a cancelation between the quadrupole associated to the anisotropic spacetime 
background and the inflation-produced one, so to explain the "quadrupole problem", occurs coincidentally (i.e., 
whatever is the direction of the axis of symmetry). The lower panel of Fig. 3 shows -P(edec) as a function of edec, and 
indicates that such a probability is not negligible for edec in the range [0.25, 1.151 x 10~ 2 . In particular, P(edcc) — 29% 
for edec — 0.65 x 10~ 2 given by Eq. (fTTj) . which amounts to have Qtot — Sobs (S. [Io|. 
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